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Abstract 

The gravitational lens equation resulting from a single (non-linear) mass concentration 
(the main lens) plus inhomogeneities of the large-scale structure is shown to be strictly 
equivalent to the single-plane gravitational lens equation without the cosmological per¬ 
turbations. The deflection potential (and, by applying the Poisson equation, also the 
mass distribution) of the equivalent single-plane lens is derived. If the main lens is de¬ 
scribed by elliptical isopotential curves plus a shear term, the equivalent single-plane 
lens will be of the same form. Due to the equivalence shown, the determination of the 
Hubble constant from time delay measurements is affected by the same mass-sheet in¬ 
variance transformation as for the single-plane lens. If the lens strength is fixed (e.g., by 
measuring the velocity dispersion of stars in the main lens), the determination of IIo is 
affected by inhomogeneous matter between us and the lens. The orientation of the mass 
distribution relative to the image positions is the same for the cosmological lens situation 
and the single-plane case. In particular this implies that cosmic shear cannot account for 
a misalignment of the observed galaxy orientation relative to the best-fitting lens model. 


1 Introduction 

Light from distant sources propagates through the matter inhomogeneities of our Uni¬ 
verse; light bundles are deflected and distorted during their propagation (e.g., Gunn 
1967a,b). This gravitational light deflection can be used to infer statistical properties of 
the large-scale matter distribution in the Universe, as was pointed out by Miralda-Escude 
(1991), Blandford et al. (1991), Kaiser (1992, 1996), Villumsen (1996), Bernardeau, van 
Waerbeke & Mellier (1996), and Jain & Seljak (1996), by measuring the mean elliptic- 
ity of the images of distant galaxies; assuming that galaxies are intrinsically randomly 
oriented, any net ellipticity is then attributed to the propagation. Estimates of the ex¬ 
pected effect have both been done numerically (e.g., Jaroszynski, Park & Paczynski 1990, 
Jaroszynski 1991, 1992, Bartelmann & Schneider 1991, Wambsganss et al. 1996) and an¬ 
alytically, using the linear or fully non-linear evolution of the power spectrum of density 
fluctuations (Seljak 1996, Jain & Seljak 1996, and references therein). 

This cosmic shear acts of course also on light bundles corresponding to gravitational 
lens systems such as multiply-imaged QSOs. Since the cosmic shear is a stochastic 
process, its value in the direction of a gravitational lens system is unknown and has to be 
determined from lens modeling. By that - as we shall see below - the cosmic shear cannot 
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be distinguished from locally generated shear, i.e., by a group of galaxies or a cluster 
close to the main lens. In addition to the shear, the large-scale matter inhomogeneities 
can also produce a convergence (which can have either sign). This convergence which is 
indistinguishable from a convergence caused by the local environment of the main lens, 
affects the determination of the Hubble constant (Kayser & Refsdal 1983, Gorenstein, 
Falco & Shapiro 1988). The primary effect here is that matter inhomogeneities between 
us and the lens change the angular-diameter distance to the lens, which essentially is the 
quantity determined from a measurement of the time delay (Narayan 1991). 

Surpi, Harari & Frieman (1996) and Bar-Kana (1996, hereafter B-K96) have investi¬ 
gated the effect of large-scale structure on multiply imaged sources. In B-K96, the effect 
of cosmic shear added to a given gravitational lens was considered and shown to affect 
the observable image positions and flux ratios. However, since the mass model for the 
gravitational lens is obtained by fitting observable image properties with model predic¬ 
tions, the addition of a cosmic shear implies that the adopted lens model is modified such 
as to reproduce the observables best. Therefore if one wants to investigate the effect of 
cosmic shear on multiple image lens systems, the combined effect of adding cosmic shear 
and modifying the lens model accordingly must be studied. 

The present paper considers some aspects of B-K96 in more detail. The main results 
to be shown are as follows: (1) B-K96 derives a lens equation including cosmic shear 
which is formally identical to the (multiple-plane) generalized quadrupole gravitational 
lens equation (Kovner 1987; Schneider, Ehlers & Falco 1992; hereafter SEF, Chap. 9). 
B-K96 shows that the analogue of the ‘telescope matrix’ (Kovner 1987) is symmetric 
if the effects of the large-scale matter inhomogeneities are considered to first order. It 
will be shown below that this matrix is manifestly symmetric to all orders. This is 
not unimportant: whereas the rms value of the cosmic shear is below 10%, it obeys 
non-Gaussian statistics, and considerably higher values are probably not rare. (2) The 
symmetry of this matrix is then used to show that the cosmic lens equation is fully 
equivalent to a single-plane gravitational lens equation, and the equivalent single-plane 
matter distribution (or, equivalently, the single-plane deflection potential) is obtained. 
(3) If a family of lens models is considered where the main lens consists of a linear 
term (a shear matrix) plus a main lens with elliptical potential, then the equivalent 
single-plane lens is contained in this family of lens models. Note that this kind of lens 
models has been applied to observed lens systems, and it appears that an (external) 
shear in addition to the ellipticity of the lens galaxy is needed in most cases where 
enough observational constraints are available (e.g., four-image systems) - see Keeton, 
Kochanek & Seljak (1996), Witt & Mao (1997). (4) The observed image positions relative 
to the observable major axis of the matter distribution in the main lens is the same in 
the cosmic lens situation and in the equivalent single-plane situation, in contrast to the 
possibility indicated by Blandford & Kundic (1996) who suspected possible misalignment 
of the images with respect to the observed orientation of the galaxy. 
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2 Light propagation in an inhomogeneous Universe 

In this section we closely follow the paper by Seitz, Schneider & Ehlers (1994; hereafter 
SSE); the reader is referred to this paper for details. 


2.1 General propagation equations 

Consider a fiducial light ray ending at the observer at redshift z — 0. The corresponding 
null geodesic is denoted by 7 q(A), where A is an affine parameter along the ray. We 
choose A such that A = 0 at the observer, and that locally A coincides with the proper 
distance. Consider a neighbouring light ray y M (A; 0) which at the observer propagates 
at an angle 0 relative to the central ray. The transverse component of the separation 
vector £ M (A; 0) — y M (A;0) — 7 q(A) is (essentially) a two-dimensional vector which shall 
be denoted by £(A; 0) (for an exact definition, see SSE). Provided 0 is sufficiently small, 
£(A; 0) will depend linearly on 0; it satisfies the Jacobi differential equation 

| // (A;0) = T(A)£(A;0) , (2.1) 

where T(A) is the optical tidal matrix evaluated at the position 7 q (A) of the fiducial ray 
at the affine parameter A; see SSE for the general definition of the optical tidal matrix. 
A prime denotes differentiation with respect to A. Equivalently, we can write 

£(A;0)=2?(A )0 , (2.2) 

which after insertion into (2.1) yields an equation for the matrix T>, 

V"(X) — T(A) V{\) . (2.3) 


From the requirement that A agrees locally with the proper distance, the initial conditions 
for £ are £(0) = 0, £'(0) = 0, or equivalently, X>(0) = 0, X>'(0) = X, where X denotes the 
unit matrix. 

Suppose the matter distribution in the Universe is characterized by a homogeneous 
matter density p(z) — (1 + z) 3 f2 0 p cr , with critical density p cr = 3H'q/( 87rG)] in that case, 
the optical tidal matrix X is proportional to the unit matrix and is given by 


r ( A) = -2(lh) a 0 (i + zfi . 

The relation between the affine parameter and the redshift z reads 


Ho (1 + z) 3 a/1 + zf2 0 ~ Ha [1 - (1 + z) 2 


(2.4) 


(2.5) 


where fl a = A/ (3Hq) is the density parameter associated with the cosmological constant 
A. In this case of a homogeneous Universe, the matrix T>{ A) remains proportional to the 
unit matrix, V( A) = D( A)X, which implies an isotropic propagation of light bundles; 
i.e., the cross-section of initially circular light bundles remains circular. The angular- 
diameter-distance D as a function of A or z can then be obtained from the differential 
equation (2.3) for the V, using (2.4). 

Adding density perturbations 5p to the Universe, the optical tidal matrix changes 
to 
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Ho 

c 


2 


( 1+^) 2 


( 2 . 6 ) 




i?o(l + ^) 5 5ij — 


(2 $,ij + 5y^,33) 


where the local coordinates are such that the light ray propagates in the ^ 3 -direction. 
The gravitational potential <I> is related to the density perturbations via the Poisson 
equation V 2( I> — AnG5p (where the differential operator is taken with respect to local 
proper coordinates) which can be justified provided the density perturbations have typical 
scales much smaller than the radius of curvature of the Universe, i.e., if that scale is 
much less than the Hubble radius c/Hq. In general, the matrix 'D(X) now attains shear 
components, and initially circular light bundles obtain elliptical cross sections. 

Finally, we add to the density perturbations a non-linear gravitational lens at redshift 
Zd, described by its surface mass density U(£j. where £ is the proper distance vector in 
the lens plane, i.e., a plane perpendicular to the light rays under consideration. If the 
fiducial light ray is chosen such that it traverses the lens plane at $, — 0 , then the 
separation vector £(A) becomes 


wn f T>(\)0 for A < A d 

^ } \ V(X)0 - V(X) [a(£(A d )) - a(0)] for A > A d 

where a(£) is the deflection angle a light ray undergoes at position £, 


(2.7) 


“(« — [ d 2 ^' £(«') ; (2.8) 

the fact that £(A) satisfies the propagation equation ( 2 . 1 ) for all A 7 ^ A d with the optical 
tidal matrix (2.6) implies that T>( A) also satisfies the differential equation (2.3). Eq.(2.7) 
is derived in Sect. 4.4 of SSE, where it is also shown that the requirement that the local 
change of direction of the light ray in the lens plane agrees with the deflection angle 
implies that the initial conditions for V are X>(A d ) = 0, X>'(A d ) = (1 + z d )Z. 


2.2 The cosmological lens equation 

In the case of an isolated mass concentration acting as a gravitational lens, there always 
exists at least one point where the deflection angle vanishes; this can be easily shown 
with the Poincare’s index theorem (see, e.g., Sect. 5.4.1 of SEF). We shall assume that 
the origin in the lens plane - and thus the fiducial ray 7 q - is chosen such that ck(0) = 0 . 

Consider sources at redshift z s , corresponding to the affine parameter A s . Let 77 = 
£(A S ), and, as before, £ = £(A d ). The lens equation in an inhomogeneous Universe then 
becomes 

7/ = Z(A s )D- 1 (A d )C-T’(A s )o : (0 , (2.9) 

where we have assumed that the lensing effect by the perturbations 5p is sufficiently small 
so that no caustic points are caused by these perturbations themselves. This implies that 
the matrices T>( A) and T>( A) are nowhere singular and thus can be inverted. In (2.9), we 
have written £ = 'D(X,])0 as independent variable; the non-singularity requirement for T> 
implies that there is a one-to-one relation between £ and 6. Eq. (2.9) maps a vector £ from 
the lens plane into a vector 77 in the source plane, just as in standard gravitational lens 
theory; the difference between (2.9) and the standard lens equation is that in the latter 
case, the matrices T> and T> are proportional to the unit matrix, i.e., they effectively 


4 



become scalars. These scalars are the angular diameter-distances, which are uniquely 
defined in terms of redshifts in a homogeneous Universe. In an inhomogeneous Universe, 
the matrices T> and T> attain shear components, so that the angular diameter-distances 
are no longer isotropic; moreover, since the density perturbations 5p form a random 
held, the optical tidal matrix along the direction to a gravitational lens has random 
components, and so there is basically no hope to obtain enough knowledge about this 
matter distribution to determine T>( Ad), T>( A s ), and T>(X a ) for a given lens system. For 
notational convenience, we define 

V d =V(X d ) ; V S = V(X S ) ; V s = t>(X s ) . (2.10) 

Furthermore it should be noted that these matrices are in general not symmetric, 
so that the cosmic perturbations induce a rotational component to the matrices T> and 
T> (this has been called ‘twist’ in SSE). However, this twist is unobservable: let (3 = 
V~ 1 (X s )r]-, then, 

0 = V- 1 V.V^(-a{()=:C ( (-a(() . (2.11) 

As we shall show in the Appendix, the matrix multiplying £ in ( 2 . 11 ) is symmetric, 
so that the matrix d/3/d$, is symmetric. The symmetry of is in complete analogy 
to the symmetry of the ‘telescope matrix’ introduced by Kovner (1987) in the frame 
of a multiple deflection gravitational lens system with at most one non-linear deflector; 
see also SEF, Sect. 9.3. The explicit algebraic proof for the symmetry of the ‘telescope 
matrix’ has been given in Seitz & Schneider (1994). We also note in passing that the 
magnification theorem (Schneider 1984; Seitz & Schneider 1992) is still satisfied: at least 
one of the images formed by the gravitational lens of any source is brighter than the 
source would appear in the same direction (i.e., with the same optical tidal matrix) in 
the absence of the lens, due to the 11011 -negativity of A7(£). However, it should be noted 
that the source can appear fainter than the same source at the same redshift would 
appear in the corresponding homogeneous Universe, since 5p can be negative; indeed, 
(Sp) — 0 by definition. 

We write the lens equation (2.11) in terms of the angle 0, by multiplying (2.11) from 
the left by Vj. and by defining rj' = T>J/3: 

ri = VjV^V s e - Dja (V d 0) =: C e 0 ~ Vjcx (V d 0) . (2.12) 

Since Cq = T>JC^D d , the symmetry of C\ implies that Cq is also symmetric. 


3 Relation to observables 

The ‘source position’ rj' is obtained from the ‘true’ source position 77 through a linear 
transformation. However, this transformation does not affect the observables. Let 6 l be 
the observed image positions, 1 < i < N, then the lens equation predicts that they have 
to satisfy 

(0 l - £F) - Vj [a, (V d 9 l ) - a (D d 0 J )] , (3.1) 

for i ^ j. a relation not containing rj'. Secondly, whereas the magnification of the images 
depends 011 the linear transformation in the source plane, the observable magnification 
ratios are unaffected by it. What can be observed are the relative magnification matrices 
Aij between images, 
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Aij — 


dr)'(O ') \ ( drj'(d r 


-l 


de 


oe 


(3.2) 


which are unchanged by a linear transformation in the source plane. Note that 


dr)' 

~d6 


C e -VjU(V d O)V d 


(3.3) 


where 

U(t) = . ( 3 . 4 ) 

Since U is symmetric, so is dr)'/86. 

We next consider the time delay function T. For that we first note that the potential 
part of the time delay is (Cooke & Kantowski 1975) 


cT po t = (1 + 2d)^(0 , (3-5) 

where 

*«) = ^rj R2 d2 «' £<«') hi ’ (3 ' 6) 

is the deflection potential, and is an arbitrary length scale; changing changes W 
by an additive constant which does therefore not affect the measureable time delays. 
From Fermat’s principle (Schneider 1985, Kovner 1990) one knows that the lens equation 
is equivalent to VT = 0, which determines T up to an affine transformation. The 
multiplicative constant is determined by that of T pot , so that 

cT = (1 + z d ) [6 ■ (C e 0 )/2 -r)'-6- V(V d 6)] , (3.7) 

and noting that d^/dO = Vja(V d 6), one sees that dT/dO = 0 is equivalent to (2.12). 
The time delay between any pair of images is then At = T ( 6 — T (0 3 ). Whereas T as 
written in the form (3.7) contains the source position, one notes that rj' is substituted 
in terms of the observed image positions using the lens equation (2.12). Therefore, the 
linear transformation in the source plane does not affect the calculation of the time delay. 


4 The equivalent single-plane lens 


We shall now show that the lens equation (2.12) is equivalent to a single-plane lens 
situation without ‘cosmic shear’. The single-plane lens equation reads 

f] = D s 6 - D da a.{D d 6) , (4.1) 


where D d , D dH . and D s are the angular diameter distances in a homogeneous Universe, 
as described in Sect. 2. We characterize quantities in the single-plane lens with a tilde, 
to distinguish them from those of the cosmological lens equation. As before, a(£) can 
be obtained from a deflection potential 


&(0 


d| 


(4.2) 
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Theorem: For every mass distribution I7(£) and matrices 'D d . T> s and t> 8 giving rise to 
the cosmological lens equation (2.12), there exists a mass distribution U(£) for which the 
single-plane lens equation (4.1) yields the same observables. 

This will be shown by construction. To simplify notation we define the matrices 


TZd — 


Vd 

D d 


V s 

K s = — 

A 


n ds 


A 

-As 


(4.3) 


which describe the deviation of the light propagation in an inhomogeneous Universe from 
that in the homogeneous one; i.e., in a homogeneous Universe these three matrices would 
become the unit matrix. Define the deflection potential 

*(€) := 1 (nd) + p ■ (b|) . (4.4) 

where B is the matrix 

B = ^ {njndn, - X) ; (4.5) 

A 

then the deflection angle becomes, according to (4.2), 

&(l) = n^a(n d l) + Bl . (4.6) 

This is now inserted into the lens equation (4.1) to yield 

t) = D S Q - D ds nJa{V d O) + D d BO . (4.7) 


If we now define r\ = (D d /D ds )r], then the lens equation (4.7) becomes after multiplica¬ 
tion with ( D d /D ds ) 

fi' = C e O-Vja(V d 0) , (4.8) 

which is obviously the same equation as (2.12), except for a linear transformation of 
the source coordinates. In particular, this single-plane lens yields the same relations 
(3.1), (3.2) and (3.7) between the image positions, the relative magnification matrices 
and the time delay function as the cosmological lens mapping. The surface mass density 
corresponding to the equivalent single-plane lens is obtained from (4.4) via Poisson’s 
equation. The first term in (4.4) then yields a contribution to I7(£) which is obtained from 
U(£) in a non-trivial way, although the corresponding potential contributions are related 
by a simple linear transformation. The second term in (4.4) yields a uniform surface 
mass density. Note that both contributions are not guaranteed to be non-negative. 

If we use the angular source position (3 = f]/D S: the lens equation takes the more 
familiar form 


0 = njndn.e 


Dds 

aT 


nJa(V d O) 


re- 


Dds 

aT 


Kja(V d O) 


(4.9) 


so that in fact the equivalent single-plane lens consists of the shear matrix T plus a term 
corresponding to the ‘main lens’. Since T is proportional to Co it is also symmetric. 

In general, the trace of T is different from 2, which means that the shear matrix con¬ 
tains an equivalent surface mass density term (which is not necessarily positive). As was 
pointed out by Gorensteiu, Falco & Shapiro (1988), adding a homogeneous surface mass 
density to a lens and at the same time rescaling the lens strength leaves all observables 
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invariant except the time delay. This mass-sheet degeneracy is therefore also present in 
the single-plane lens and is not a novel feature of the cosmological lens equation. 


5 Special case: Elliptical isopotentials 

We now consider the case that the lens is described by elliptical isopotential curves, so 
that 

V(S) = F(u) , (5.1) 

with 

u = £ 2 (l - ecos[2 (<p -#)]) ; (5.2) 

here, £ has been expressed in polar coordinates £ = £(cos <p, simp), e describes the ellip- 
ticity of the potential, and d is the direction of the major axis of the potential. This kind 
of lens model has been introduced by Blandford & Kochanek (1987) and investigated 
in detail by several authors (e.g., Kassiola & Kovner 1993). The corresponding mass 
distribution has no elliptical isodensity curves, and depending on the shape of the radial 
profile F(u) and the ellipticity e, dumb bell shaped isodensity curves or even negative 
surface mass densities are obtained from . Leaving aside these potential difficulties, 
models of the form (5.1) have often been used for modeling gravitational lens systems. 
Whereas the surface mass density is not elliptical, one can still define a major and minor 
axis, which coincide with the major and minor axes of X V. 

Expressed in terms of the observable 0, W(T) ( \6). the observed major axis of the lens 
does not coincide with 0, but is rotated by the matrix 'D ( \. But the same is true for 
the equivalent single-plane lens, see (4.4). Note that the first term in (4.4) is again an 
elliptical potential in however with different ellipticity and orientation, and a modified 
radial profile. However, the matrix which rotates \P(DaO) with respect to 'T is the same 
as that rotating the intrinsic lens direction 0 into the observed one. Hence, the observer 
in the inhomogeneous Universe sees the same orientation of the lens as the corresponding 
observer in the homogeneous Universe sees of the equivalent single-plane lens. Since the 
image positions G l are the same in both cases, this implies that the positions of the images 
relative to the observed orientation of the lens are the same in both cases, in contrast 
to the conjecture made in Blandford & Kundic (1996). In fact, not only is the observed 
orientation the same in both situations, but the potentials do agree in both cases when 
expressed in terms of G\ Note that this is not the case for the surface mass densities. 

Finally, it can be easily seen that if one considers a family of models where the 
primary lens has elliptical potential curves and in addition a shear matrix is included, 
then the cosmological contributions Cg and 77<i yield an equivalent single-plane lens which 
is included within this family. 


6 Summary and discussion 

In this paper the effects of large-scale matter inhomogeneities on gravitational lensing by 
galaxies were considered. The lens equation in an inhomogeneous Universe was derived; 
the only approximation entering (2.9) is that the length scale of the matter inhomo¬ 
geneities is much larger than the linear size of a light bundle which encloses all light rays 



from the source to the observer (i.e., for a typical multiply imaged QSO, this size is of 
order 20kpc). In that case, light propagating from the source to the lens and from the 
lens to the observer can be described by the Jacobi differential equation. Note that for 
the validity of (2.9) one does not need to assume that the inhomogeneities are weak - 
or in other words, that the matrices 77a, 77 s and 77,i s are close to the unit matrix. The 
analogy of the ‘telescope matrix’ which was defined by Kovner (1987) has been shown 
to be symmetric in general. This symmetry property is essential for showing that for 
arbitrary propagation matrices 77,i, 77 s and 77<i s caused by the large-scale structure, there 
exists an equivalent single-plane gravitational lens such that all observables are the same 
- image positions, relative magnification matrices and time-delays. However, in general 
this equivalent single-plane lens contains a homogeneous matter sheet, not necessarily 
with positive surface mass density. Therefore, the mass-sheet degeneracy as discussed in 
Gorenstein, Falco & Shapiro (1988) does apply. 

For a single lens plane, this mass-sheet degeneracy can be broken if independent 
observational information about the lens can be obtained. For example, if the lens is well 
described by an isothermal mass distribution, the measurement of the velocity dispersion 
determines the lens strength and thus fixes the mass-sheet transformation. If this is done, 
the time delay then depends on the matrix 77a, and thus matter inhomogeneities between 
us and the lens affect the determination of the Hubble constant. This fact was noted 
by Narayan (1991) who showed that the measurement of the time delay determines the 
angular diameter distance to the lens, which is affected by 77a, in agreement with B-K96, 
where the size of this effect has been calculated using the non-linear evolution of the 
power spectrum. Thus, lensing by the large-scale structure does affect the determination 
of the Hubble constant, although the magnitude of the effect (i.e., the rms deviation of 
77a from the unit matrix) is of the order of a few percent; see Fig. 2 of B-K96. Therefore, 
in agreement with the conclusion of Surpi, Harari & Frieman (1996), cosmic shear does 
not seriously compromise the lensing method for the determination of Hq. 

Whereas the mass distribution of the equivalent single-plane lens depends non- 
trivially on the lens mass distribution and the propagation matrices 77, it has been shown 
that if the lens is chosen to have elliptical isopotentials then the equivalent single-plane 
lens consists of a main component with elliptical isopotentials plus a shear matrix. Fur¬ 
thermore, if a family of lens models is considered, consisting of an elliptical isopotential 
component and a shear matrix, then the equivalent single-plane lens is also within this 
family. For this case it was shown that the image positions relative to the orientation 
of the mass distribution as seen by an observer is the same in the inhomogeneous Uni¬ 
verse and in the single-plane lensing situation. Note that lens models consisting of an 
elliptical lens (either elliptical isodensity contour or elliptical isopotentials) plus external 
shear seems to be required for those lens systems in which there is a large number of 
observational constraints (e.g., Keeton, Kochanek & Seljak 1996, Witt & Mao 1997). 

I thank S. White, J. Frieman, and in particular S. Seitz for valuable discussions, 
and M. Bartelmann for comments on the manuscript. This work was supported by the 
“Sonderforschungsbereich 375-95 fur Astro-Teilchenphysik” der Deutschen Forschungs- 
gemeinschaft. 
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Appendix 

We shall show in this Appendix that the matrix multiplying in (2.11) is symmetric. 
To do this, we define the matrix 


A(A) :-D(A)D- 1 (A d ) 


(Al) 


and the matrix 

C'(A) - V-\X)V(X)V-\X d ) = V~\X)A(X) , (A2) 

in both cases for A > A d . The proof proceeds in the following four steps: 

(1) The matrix A(A) satisfies the differential equation (2.3), which implies that A 
and D satisfy the same differential equation. This can be trivially shown by inserting 
(Al) into (2.3). 

(2) For any two solutions X and Y of the differential equation (2.3), one finds 

-i[(X T ) / y-X T W] =0 , (A3) 

as can be easily verified by straight insertion into (2.3). This relation implies that 

5(A) :=P'(A)P“ 1 (A) (A4) 


is symmetric. This can be shown as follows: Eq. (A3) implies, by setting X — Y — that 
(TY-yV — fYf)' is constant. Due to the initial condition at A = A d , X> T (A d ) = 0 = X>(A d ), 
this constant is the zero matrix, so that ( V l )'T> — T> l T>' for all A > A d . In particular 
this implies 

(V T y = V^V'V- 1 . (A5) 

To show the symmetry of 5, we consider 


S — S T — v'v 


-1 


V 


-1 


(v T y = v’v 


-i 


V 


-1 


V l V'V~ l = 0 


(A6) 


where we have used (A5) in the second step, and in the final step the fact was employed 
that the inverse and transpose operations do commute. Hence, S is symmetric. 

(3) If we set X — A and Y = T>, we find that the constant in the bracket of (A3) is 

(A t )'P- A t V' = ~(l + z d )l , (A7) 

due to the intial conditions A(A d ) = X, X>(A d ) = 0, X>'(A d ) = (1 + z d )X. Hence, 

(1 + z d )V~ l - A t S - (A T )' . (A8) 


(4) The final step is taken by noting that liniA^Ad [(^ — -^d)C < (A)] is symmetric. We 
now show that C' is symmetric for all A > A d , which then proves that C is symmetric 
for all A — A d , and in particular this is true for = C(A S ). Using 

VV~ l =X ^ (V- 1 )'=-T)- 1 S , 


as follows from differentiation, we find in turn: 
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C' = V~ X A! + {V~ l )'A = V^A' - V-'SA 

1 


(. A t SA ' - ( A t )’A' - A t SSA + (A t )'SA) 


m 


(1 + z<±) 

However, this final expression is manifestly symmetric, which completes the proof. 
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